SCIPP 93/21 
July 1993 



Phenomenology of Gluino Searches at the Tevatron 



Howard E. Haber 

Santa Cruz Institute for Particle Physics 
University of California, Santa Cruz, CA 95064 

Abstract 

Present data indicates that the gluino (if it exists) must be heavier than about 
95 GeV. During the next few years as the Tevatron integrated luminosity increases, 
gluino searches will be able to probe the mass range between 100 and 200 GeV. For 
masses in this range, a variety of gluino decay modes can provide viable signatures 
for gluino detection. Apart from the classic missing transverse energy signal, the 
detection of high transverse momentum like-sign dileptons may be the cleanest 
signature for gluino production. Other signatures such as the production of a hard 
photon in the gluino cascade decay may also play an important role in confirming 
the supersymmetric origin of events originating from gluino production and decay. 
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1. The Gluino Mass Limit, circa 1992. 

The Tevatron has recently completed its 1992-93 run and has delivered about 
25 pb" 1 of data to the CDF and DO detectors. Next year, the Tevatron hopes 
to triple the accumulated luminosity of the 1992-93 run. With this large increase 
of data from the highest energy collider in the world, there is some hope that 
evidence for deviation from the Standard Model could soon be at hand. At this 
meeting on the status of supersymmetry in 1993, it is appropriate to reflect on the 
prospects for the discovery of supersymmetry at the Tevatron in the near future. 
Since squarks and gluinos have the largest production cross-sections at a hadron 
collider, it is these particles that have attracted the major attention of the two 
detector Collaborations in their searches for supersymmetry. 

In supersymmetry searches at hadron colliders, two basic scenarios emerge. 
Gluinos and squarks are produced in pairs, so that possible final states are gg, 
qq and gq. If Mq < Mg, then g — > qq (or qq), and one should first concentrate 
on the signatures of squark production at the Tevatron. Likewise, if Mg < Mq, 
then q — > qg, and it is appropriate to first focus attention on gluino production. 
Supersymmetric models do not yield a definitive statement as to which case is more 
likely to be realized in nature, although there may be a slight theoretical bias in 
favor of the lighter gluino. In this paper, I choose to consider the case of Mg < Mq 
and examine the phenomenology of gluino production and decay at the Tevatron. 

During the past decade, many experiments have conducted searches for 
gluinos. 1 The results of all such searches have been negative. From these nega- 
tive searches, gluino mass limits have been obtained. In the literature, there has 
been some disagreement on the extracted gluino mass limits. A comprehensive 
discussion of these disagreements lies beyond the scope of this paper. Instead, I 
will state my personal opinions, and indicate how I draw my conclusions. A brief 
summary on the status of gluino searches is presented in Table 1. 

Two main questions have been debated: (i) are very light gluinos allowed? and 
(ii) what is the maximum gluino mass ruled out by collider data? It is convenient 
to separate the discussion into two mass regimes: Mg < 3 GeV and Mg > 3 GeV. 
The borderline between these two regions has been called the "light gluino window" 
in the literature. Recently, there has been a revival of interest in this region. A 
light gluino has been advocated in order to explain a possible inconsistency in the 
determination of a s by different techniques. 9 Of course, experimental searches will 
be the final arbiter of the existence or non-existence of the light gluino window. 

A variety of beam dump experiments and beam contamination searches rule 
out gluinos with Mg < 3 GeV, with some dependence on the gluino lifetime. 2 



| Here, I do not distinguish between squarks and their antiparticles. On the other hand, gluinos are 
self-conjugate, a property that plays an important role in section 4. 
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Table 1. Gluino Mass Limits 



M~ Sensitivity 


Experiment 


Comments 


M~ < 1-10 GeV 

9 


A variety of beam dump 
experiments and beam 
contamination searches 2 


Neutron backgrounds are prob- 
lematical if Mg ~ m n ; mass 
limits depend on gluino lifetime 


M~ & 1 GeV 

9 


ip -> IVg 


rig — ^ 53 bound state; 


M~<,3 GeV 




color factor of 27/4 enhance- 
ment over production of r)Q 


3 <> M~ <, 4 GeV 


Based on UA1 data 


Analysis of ref. 4; important use 
of ggg production 


4 <, M~ <, 53 GeV 

9 


UA1 5 


Assumes g — > qqj with 


16 ^ M~ <, 79 GeV 


UA2 6 


My £ 20 GeV 


30 ^ M~ <, 95 GeV 


CDF 7 


Upper limit includes effects of 
gluino cascade decays with 
,u = -250 GeV and tan j3 = 2 


M~ < 132 GeV 

9 


ref. 8 


Based on CDF gluino mass lim- 
its and LEP limits on chargino 
and neutralino masses [assuming 
the validity of eq. (2.1)]. 



Perhaps the most definitive negative result is the non-existence of a pseudoscalar 
gg bound state (sometimes called the rjg). Such a state, if it existed, would be 
prominent in radiative quarkonium decay. 10 Because gluinos are color octets, the 
decay rates for ip — * 777^ and T — > 7?^ are enhanced by a color factor of 27/4 
over the radiative decay into 7^q, where t]q is a pseudoscalar meson bound state 
of color-triplet quarks (e.g., r/, rj', and r) c ). No such prominent pseudoscalar state 
is seen in the data. The published limits of the CUSB collaboration 11 at CESR 
rules out gluino masses between 0.6 GeV and 2.2 GeV. Subsequent analysis of the 
CUSB data 3 extended the upper mass limit to 2.6 GeV. Presumably, additional T 
decay data now exists that could push the gluino mass bound above 3 GeV. 

For gluino masses above 3 GeV, one must turn to the collider data. Gluino mass 
limits have been presented by the UAl 5 and UA2 6 Collaborations based on CERN 
pp Collider data and by the CDF Collaboration 7 based on Tevatron data. These 
analyses assume that the gluino decays inside the detector with the emission of the 
lightest supersymmetric particle (LSP) which carries off undetected ("missing") 
transverse energy. Only the UAl data is sensitive to light gluino masses. The 
UAl limits rule out gluinos in the range 4 < Mg < 53 GeV. This apparently 
leaves open the light gluino window. However, I believe that the UAl limits are 
too conservative. In particular, the UAl analysis did not take into account the 
mechanism gg — > ggg which is especially important for Mg < 5 GeV. In particular, 
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the kinematical configuration in which a hard gluon recoils against the two gluinos 
which are emitted in the same hemisphere can lead to substantial missing energy 
and pass the UAl cuts. Barnett, Kane and I investigated this mechanism carefully 4 
and concluded that the UAl data definitively ruled out gluino masses as light as 
3 GeV. This essentially eliminates the possibility of the light gluino window. 

The analyses of the UA2 and CDF Collaborations are not sensitive to the light 
gluino. However, their excluded mass regions overlap those of the UAl analysis and 
therefore extend the UAl gluino mass bound. The UA2 Collaboration concluded 
that Mg < 79 GeV. The CDF Collaboration has extended the gluino mass bound 
further; published results 7 have been presented based on 4.3 pb^ 1 of data from the 
1989-91 Tevatron runs. In order to properly interpret the CDF data and extract 
gluino mass limits, one must incorporate all the allowed gluino decay patterns in 
the analysis. Typically, one assumes that the LSP is the lightest neutralino (xi)- 
Then, the simplest gluino decay is g — » qqxv sometimes called "the direct decay 
to the LSP". This is expected to be the dominant decay for light gluinos* For 
heavier gluinos, decays into heavier neutralinos and charginos 

(1.1) 

9 -»■ QQ Xi 

become kinematically allowed, in which case, the branching ratio for the direct 
decay to the LSP drops significantly below 1. This is relevant for phenomenology 
since it is the direct decay into the LSP that produces the largest missing energy 
signature. When the gluino decays first into heavier neutralinos or charginos, these 
particles subsequently decay, eventually producing the LSP at the end of the decay 
chain. The end result is a gluino cascade decay that produces less missing energy 
than the direct decay to the LSP. 

The CDF Collaboration has presented 7 gluino mass limits under the assump- 
tion that BR(g — * qqxi) — 1- However, such mass limits are not very meaningful, 
since they correspond to a very unlikely choice of parameters for the supersym- 
metric model. For a more realistic example, consider the minimal supersymmetric 
extension of the Standard Model (MSSM), with gaugino Majorana mass parame- 
ters related according to the unification relation 

Mi/Mj = a-ijoLj , (1.2) 

where the % labels the gauge group, ctj = gf/47r, and g\ = (5/3)</ 2 . In this case, 

* There is an implicit assumption here that the gluino is heavier than Xi- Otherwise, the gluino 
would be the LSP and hence stable. Note that the CUSB limits quoted above also apply to stable 
gluinos. Beam contamination experiments then rule out stable gluinos for Mg < 10 GeV. 2 Heavier 
stable gluinos are disfavored for both theoretical and cosmological reasons, although this possibility 
cannot be completely excluded. 



4 



the neutralino and chargino masses and mixing angles are determined by three 
parameters: the gluino mass (Mg), the ratio of Higgs vacuum expectation values 
(tan/?) and the supersymmetric Higgs mass parameter (/i). 12 Thus, for a given 
value of gluino mass, the parameters /i and tan (3 determine the relative branching 
ratios of the gluino into charginos and neutralinos [eq. (1.1)]. For gluino masses 
now being probed at the Tevatron, the probability of the cascade decays (in which 
the LSP is produced at the end of a multistep decay chain) is typically larger 
than the probability of direct decay to the LSP. Thus the missing energy signal is 
degraded and the true CDF gluino mass limit should be less than the quoted limit 
for BR(g — > qq^{) — 1. In ref. 7, the CDF Collaboration has included the effect 
of the cascade decays on their analysis. For one typical set of parameters, they 
deduced an upper gluino mass bound of about 95 GeV. This bound will surely be 
improved (by both the CDF and DO Collaborations) once the 1992-93 Tevatron 
data is analyzed. 

The bound Mg > 95 GeV is the best bound presently available based on direct 
gluino searches. However, under the assumption of the unification of gaugino 
mass parameters [eq. (1.2)], the search for neutralinos and charginos at LEP can 
indirectly lead to a bound on Mg. For example, in ref. 8, Hidaka concludes that 
Mg > 132 GeV, based on the LEP limits on neutralino and chargino masses and 
including the effects of cascade decays on the CDF gluino search analysis mentioned 
above.^ 

To conclude this introduction, my personal opinion is that based on current 
published data, Mg ^125 GeV, with no open light gluino window. Looking to the 
near future, the Tevatron runs in 1992-93 and 1993-94 will yield about 25 pb _1 
and 75 pb _1 per detector, and should increase the gluino mass sensitivity to about 
200 GeV. The purpose of this talk is to indicate some of the phenomenological 
methods available to achieve this sensitivity. In section 2, I present a brief review 
of the theory of gluino cascade decays. For gluino masses that can be detected 
at the Tevatron, there are three important gluino decays: (i) the direct decay to 
the LSP (g — > (ii) g — > qqx® > anc ^ U ) 9 ~^ QQ'Xi- m case (")> a viable 

signature results if x\ — > ^7; this is the subject of section 3. In case (iii), note 

that the Majorana nature of the gluino"'' leads to equal probability for producing 
either sign chargino. Thus, if the chargino decays semi-leptonically, gg events can 
lead to like-sign dileptons in the final state. This is the subject of section 4. 



f The inferred gluino mass limit of ref. 8 was based in part on a preliminary CDF analysis of their 
1989-91 data, and should probably be reduced slightly in light of the published CDF results. Note 
that the LEP limits on chargino and neutralino masses are basically at their kinematic limits and 
hence will not change with further running at y/s = mz- 

t Here, I am using a broader definition of Majorana to include neutral particles that transform under 
real representations of the underlying Standard Model gauge group. 
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2. The Theory of Gluino Cascade Decays 

At the Tevatron, the dominant production mechanism for gluinos is gg — > gg 
via s-channel gluon and t- and -u-channel gluino exchange. 13 The Iggg coupling 
follows from QCD gauge invariance, so that the gluino cross-section depends on 
only one unknown parameter — the gluino mass. Gluinos can also be produced in 
qq annihilation via squark-exchange. However, this produces only a minor change 
in the overall gluino cross-section. Thus, gluino phenomenology depends on the 
gluino mass and the parameters that govern the gluino decay branching ratios. 
The most important gluino decay modes are tree-level three body decays into a 
neutralino or chargino and a qq pair [eq. (1.1)]. With the exception of g — > qc${, 
the neutralino or chargino in the final state will decay into a lighter neutralino or 
chargino and so on until the LSP (assumed to be x^) is produced. Branching ratios 
for the gluino cascade decays depend on various supersymmetric parameters. For 
convenience, I list below the parameters of the MSSM: 12 

1. The gaugino Majorana mass parameters (Mi, M<i and Ms). 

2. The supersymmetric Higgs mass parameter, \i. 

3. The ratio of Higgs vacuum expectation values, tan/5. 

4. The mass of the CP-odd Higgs scalar, m^o. 

5. Diagonal (fi-fl and Jr-Jr) squark and slepton mass parameters. 

6. Off-diagonal (/l-/r) squark and slepton mass parameters (which depend on //, 
tan/3 and various U A" -parameters.) 

The gluino mass is given by Mg = \Ms\, and the chargino and neutralino masses 
and mixing angles are determined by Mi, M2, \i and tan/3. The Higgs sector 
parameters are fixed by tan/3 and m^o. This is relevant to the considerations 
here since neutralinos and charginos can decay into final states containing Higgs 
bosons. Finally, gluino production rates and branching ratios depend very weakly 
on squark masses. (Here, I assume that Mg < Mg and gq and qq production 
rates can be neglected.) As mentioned above, the gg cross-section is dominated 
by gluon and gluino exchange. The gluino decays via squark exchange, but this 
just means that the squark mass determines the overall normalization of the gluino 
decay width. To the extent that squarks are roughly degenerate in mass, gluino 
branching ratios will be independent of the squark mass. For heavy gluinos relevant 
for Tevatron searches, the gluino decays essentially instantaneously (with no visible 
gap between production and subsequent decay). In considering the subsequent 



The coupling of gluinos to qq is independent of the flavor of q. Moreover, with the exception of 
the top-squark, all squarks are expected to be roughly degenerate in mass. Interesting effects could 
arise if top-squark mixing effects are important. However, the decay g —> ttx\ is not kinematically 
allowed for gluinos accessible to Tevatron searches. Thus, it is probably safe to ignore the variation 
of squark masses. 
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decay of charginos and neutralinos in the gluino decay chain, all the supersymmetric 
parameters enter to some extent. For example, charginos and neutralinos produced 
in gluino decays at the Tevatron will typically decay via squark, slepton, W or Z 
exchange into three-body final states. However, even in this case, the sensitivity to 
squark and slepton parameters is not significant, since squark and slepton masses 
almost certainly lie above mz- 

Two critical assumptions underlie the analysis presented in this paper. First, 
as discussed in section 1, I shall assume that the gluino is lighter than all squark 
masses (with the possible exception of the top-squark). As noted above, once 
I assume that M g < Mg, the precise value of Mq is not important. Second, I 
shall impose the grand unification relation which states that the gaugino mass 
parameters Mi, M 2 and M3 are all equal at some very large grand unification (or 
Planck) scale. I then use one-loop renormalization group evolution to determine 
the value of the Mj at the electroweak scale. One finds that the scaling of Mj 
is proportional to the squared coupling constants q.{ [eq. (1.2)]. As a result, the 
low-energy values of M\ and M 2 can be expressed in terms of the gluino mass 

M 2 = (g 2 /g 2 s )Mg ~ 0.285M^ , Mi = (5g' 2 /3g 2 )M 2 ~ 0.483M 2 . (2.1) 

This relation underlies nearly all phenomenological analyses of the MSSM that 
appear in the literature. Perhaps it is time to begin to consider some of the 
phenomenological consequences of the violation of eq. (2.1). However, I will not 
pursue this alternative here. 

Barnett, Gunion and I have studied the branching ratios for the gluino, neu- 
tralino and chargino as a function of the supersymmetric parameters. 14-16 (See 
also the results of ref. 17.) Below, I summarize some of the highlights of our 
investigation. 

1. The gluino branching ratio for the direct decay to the LSP is small over a 
substantial portion of parameter space. This conclusion relies on the rela- 
tion among the gaugino mass parameters [eq. (2.1)]. In particular, BR(g — > 
qqXi) < 0.14 for M § £ 500 GeV. For gluino masses accessible to future Teva- 
tron searches, this branching ratio may be somewhat larger (depending on the 
values of /x and tan/3), although it is less than 50% in all but a narrow region of 
parameter space, as shown in fig. 1. Thus, most gluino decays will be cascade 
decays, which dilute the famous missing transverse energy signature. 

2. The dominant gluino decay mode [except for regions of parameter space where 
tan/5 is near 1 and ^ 0(mz)] is g — > X\ + X where X is either ud or cs 
(or the charge conjugate pair, depending on the sign of X\)- For example, if 
tan P — 4 and M g = 120 GeV, the sum of the branching ratios for gluino decay 
into xf is approximately 58% (with very weak dependence on /x; see fig. 1). 
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Moreover, because the gluino is a Majorana fermion, 

Y(g^udxt) = T(g^udx 1 ). (2.2) 

Thus, in gg production, like-sign charginos can be produced, which can lead to 
like-sign dileptons if both charginos decay semi-leptonically. This is the subject 
of section 4. 

3. The decay g — > qqx% is competitive with the direct decay into the LSP. It is the 
second most important gluino decay mode over a substantial range of MSSM 
parameter space. 

The phenomenology of gluino cascade decays 18 will depend in detail on the 
neutralino and chargino branching ratios. Since we are interested here in gluinos 
in the mass range between 100 and 200 GeV, the corresponding masses of X\ and 
X2 will be such that tree- level two body decays X\ ~^ W^X? and x\ — * Xi^ are 
kinematically forbidden. If no other tree-level two body decays are possible, then 
the dominant decays of xf and x° will be: xf -> Xiff and x° -»• Xi//, where / 
is either a quark or lepton. Although this is certainly the most likely situation, it is 
important to consider two alternative cases. If neutral Higgs bosons (the CP-even 
h° and the CP-odd A ) are sufficiently light, then the two-body decays \\ ~* 
and %2 ~~ * X?^° could be kinematically allowed. For example, ref. 16 argued that 
a significant region of parameter space exists in which x® ~^ is the dominant 
decay mode. Since that analysis, two new pieces of information suggest that the 
Higgs decay mode of xj] will not be significant at the Tevatron gluino search. First, 
the LEP Higgs search implies that m^o 60 GeV. 19 Second, radiative corrections 
can significantly increase the theoretical expectation for m^o. 20 As a result, the "3- 
body decay region" discussed in ref. 16 expands substantially at the expense of the 
"light-Higgs region". Finally, although I have assumed in this paper that squarks 
are heavier than gluinos, it could happen that sleptons are significantly lighter 
than the squarks (as suggested by renormalization group evolution of low-energy 
supergravity models 21 ). If the sleptons are lighter than the lightest chargino or 
neutralino, then two-body decays such as xf ~^ xf ~^ X2 — y ^> and/or 
%2 —> vv may be kinematically allowed. We shall see in section 4 that if these 
two-body chargino decay modes are present, then the like-sign dilepton signal in 
'gg events is significantly enhanced. 

Finally, there is one more two-body decay mode of potential significance — the 
radiative decay x\ ~ * X?7- 22 This is not a tree-level process; it occurs only at one- 
loop. Nevertheless, it turns out that the radiative decay mode can compete with 
the three-body modes (and in some cases it can be the dominant mode) in certain 
parameter regimes. This is the subject of section 3. 

We have seen above that the branching ratio for g — > qqx^i 1S expected to be 
small (less than 50%, and more likely closer to about 20%). As a result, the missing 
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transverse energy signal is degraded relative to the case of a 100% branching ratio 
for the direct decay to the LSP. Nevertheless, the missing transverse energy signa- 
ture will still be important. 23 For example, events in which one g decays directly 
to the LSP and the second g decays arbitrarily will still produce a significant miss- 
ing transverse energy signal. However, in order to maximize the efficiency of the 
Tevatron gluino search, other gluino signatures should be examined. In particular, 
we noted above that for 100 < M~ g < 200 GeV, 

BR(g - xfqq') > BR(g - ?&qq) > BR(g - XiQq) (2.3) 

over nearly all of the relevant MSSM parameter space. Thus, it is important to 
examine distinctive gluino signatures arising from the g decay to xf an d X*2- 

3. Hard Photons in Gluino Cascade Decays 

Consider gg events in which one gluino decays to x\ v ^ a 9 ~^ qq\\-, while the 
decay mode of the second gluino is arbitrary. In this section, I shall consider the 
case where \\ — > x\l- Such events are characterized by the emission of a hard 
photon, one or more large Et jets, and significant missing Et- Events of this type 
are already being searched for at the Tevatron. 24 In this section, I shall present 
the preliminary results of a calculation to estimate the fraction of 'gg events with 
a hard photon. 25 

The branching ratio for gluino decay into \\ is typically between 20% and 30% 
for gluino masses of interest to the Tevatron search. We then need to compute 
BR(x2 ~^ Xi7) as a function of the MSSM parameters. The complete amplitude 
for the one-loop process x% — > ^7 has been computed in ref. 22. To evaluate the 
branching ratio, one must also compute the sum of all tree- level decays of the x.2- 
It is clear that any kinematically allowed tree-level \\ decay into two-body final 
states would dominate the one-loop radiative decay (as well as tree-level decays 
into three-body final states). For the range of neutralino masses relevant for the 
Tevatron gluino search, the only possible allowed two-body (tree-level) decays are 
X2 —* Xih° an d X2 ~^ However, as argued in section 2, if one takes into 

account the LEP Higgs mass limits and the results of the one-loop corrected MSSM 
Higgs masses, one would conclude that the tree-level x\ decays to hP and/or A are 
likely to be present only over a small region of parameter space. Henceforth, I shall 
assume that these decays are kinematically forbidden. In this case, the three-body 
(tree-level) X2 decay modes are dominant. However, I shall now show that the 
radiative decay X2 ~^ X?7 is competitive over a rather large region of parameter 
space of interest to the Tevatron gluino search. 
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In figs. 2-4, I have plotted BRixi ~^ Xi7) as a function of MSSM parameters 
for a variety of gluino masses. These figures show that in the region of negative //, 
the typical value of BR(j& -> x?7) is 10%. This means that in the Tevatron gluino 
search, approximately 5% of all gg events should contain a hard photon originat- 
ing from %2 — > x?7. It remains to be seen whether the reducible backgrounds can 
be removed (e.g., events in which a neutral jet or leading 7r° is misidentified as a 
photon) thereby leaving a viable signal. I suspect that the gluino will not be dis- 
covered via the radiative \2 decay. However, if evidence for the gluino is uncovered 
in other channels, the detection of hard photon events could help to confirm the 
gluino interpretation of the other signals as well as pinpoint some of the MSSM 
parameters. 

4. Discovering Gluinos with Like-Sign Di-Leptons 

In this section, I shall focus on the striking experimental signature of two 
isolated leptons which can arise from gluino pair production. Half of the events of 
this type will have leptons with the same sign of electric charge. This signature, 
which is analogous to the opposite-sign lepton signal for the top quark, may yield 
sensitivity much superior to the missing-energy gluino signature at the Tevatron, 
depending on the parameters of the supersymmetric model. The results of this 
section are based on work in collaboration with R.M. Barnett and J.F. Gunion 26 
(see also ref. 27 for an earlier version of this work). Complementary work that also 
considered the like-sign dilepton signature can be found in refs. 18 and 28-33. 

As discussed in section 2, for gluino masses accessible to the Tevatron search, 
the dominant gluino decay mode is g — * qqxi over a large range of supersymmetric 
parameter space. Leptons can result from decays such as X\ ~^ W^Xi ^ ±Z/ Xi) 
where the W ± is either on-shell (if kinematics allow) or off-shell. If the sleptons 
are lighter than the chargino, then the two-body decays xf ~^ ~^ ^ ±v x\ 
and X\ ~ * vl^ — > i^v\\ are allowed and may dominate. Since the gluino is a 
Majorana fermion, it has the distinctive property of decaying with equal probability 
into fermions and antifermions. Thus, an excellent signature for pair production 
of gluinos results from events in which both gluinos decay to a chargino of the 
same sign, yielding like-sign dileptons (£ + £ + or in the final state. The 

probability for the production of like-sign and opposite-sign leptons is equal, and 
the characteristics of the two classes of final states are identical. Observation of 
this distinctive result would be extremely helpful in identifying the origin of the 
events. 



* The LEP limits on neutralino and chargino masses tend to rule out regions of —40 GcV ^ fj, ^ 
80 GeV (with some variation depending on Mi and tan 0). 19 
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To be specific, we have taken the branching ratio for gluino decay to the lightest 
chargino to be BR(g — > qq'x\) — 0-58, a result which holds to good accuracy for 
all \fx | ^ Mg/3, mz- 15 Moreover, this value for the branching ratio is approximately 

valid for nearly all MSSM parameters of relevance to the Tevatron gluino search? If 
the £ and v are heavier than the chargino, the chargino decays dominantly into the 
LSP plus a real (or virtual) which then decays 22% of the time into electrons 
and muons. Thus the branching ratio for the decay chain g — > qqt vx\ is likely 
to be as large as 13%, with equal probability to produce a lepton of either sign. 
(For simplicity, the r-lepton will be neglected from our considerations.) Since 
gluinos are produced in pairs, the number of dilepton final states resulting from 
the decay of the two gluinos would be about 1.6% of all gg events, of which half 
would have a pair of like-sign leptons. However, if Mj and/or M~ < M~±, and 

M~± < mw + then the two-body decays of the xf to tv and/or vt will 

have approximately 100% branching ratio into leptonic final states (including the 
r-lepton). If we neglect r-leptons, we find BR(g — > qql ^Xi) close to 40%. A 
remarkable 15% of all gg events would yield a dilepton final state. 

Thus, we propose that the Tevatron search for events with hadronic jets (two 
from each gluino), missing energy due to the LSP and neutrinos in the final state, 
and a dilepton pair which can come in one of the following like-sign combinations: 
e e , e fj, , fi^fi^, and the corresponding opposite-sign combinations. The events 
would be very similar to those arising from tt production in which the leptons come 
from primary decays of the t and i. Thus, distinguishing the source of opposite- 
sign events might be difficult. Because the efficiency for tagging 6-jets is low and 
because some fraction of gg events will contain 6-jets from a hard radiated gluon, 
6-jets may not be a useful tool for separating tt from gg events. For this reason 
we will focus on like-sign dilepton final states for which tt production yields a 
background only through t — > bl~v and t — > bX, b — > cl~v (or the corresponding 
charge-conjugated decay chain). This background would be quite small since the 
lepton from the 6 decay would very rarely be isolated. 

We have evaluated the rates for the dilepton signal in gluino pair production 
at the Tevatron. In order to roughly account for realistic experimental conditions, 
we have employed a parton-level Monte Carlo, which included resolution smearing 
but no fragmentation, to model the 'gg events"'" The surprisingly large potential 
for gluino discovery at the Tevatron becomes apparent by giving the number of 
dilepton (opposite- plus like-sign) events obtained in a 25 pb _1 year. For the case 



f For example, when M g Ss 200 GeV and tan/3 <; 4, BR(g — > qq'xf ) varies between about 45% and 
65% as n is varied. 

J A description of a similar program to analyze the characteristics of supersymmetric events at the 
CERN pp collider can be found in ref. 4. 
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Table 2. Number of l^i^ plus i^t^ events after lepton cuts for various M g values at the Tevatron 

with an integrated luminosity of 25 pb _1 . The various decay modes of the X\ are indicated by W 

(WXi), I {iv), v (vi). These rates assume the branching ratios quoted earlier (which are, in fact, 
typical over a wide range of parameters). 



Mode 


W 


W 


W 


I 


V 


V 


I 


I 


I 


V 


V 


V 


I 


V 


M(g) 


140 


140 


140 


160 


160 


160 


180 


180 


180 


180 


180 


180 


200 


200 


M(xf) 


80 


60 


45 


60 


80 


60 


80 


80 


45 


80 


80 


60 


80 


80 


M(£ or v) 








40 


70 


50 


75 


55 


40 


40 


60 


50 


75 


40 


Events 


9 


6 


3 


22 


9 


12 


30 


22 


9 


28 


20 


7 


15 


15 



where the X\ decays to Iv or £u, the net branching ratio of 15% quoted above 
yields roughly 1440, 476, 183, 79, 37, and 18 dilepton events (before cuts) for 
gluino masses of 100, 120, 140, 160, 180, and 200 GeV, respectively. Dilepton rates 
originating from X\ ~^ W^Xi decays are a factor of roughly 9 smaller. 

To estimate the rates for detectable dilepton events at the Tevatron, we assume 
a trigger which requires that the leading (secondary) lepton has Ej* > 15 (10) 

GeV. In addition, we require \r)\ < 2.5 and isolation for both leptons. Associated 
hadronic jets are not required. The probability that events pass these cuts de- 
pends in detail upon the masses of the particles involved in the decay chains, but 
is roughly 50% for much of parameter space. We have explored a range of super- 
symmetric parameters for which the mass M~± varies between about 45 GeV (its 

current lower bound from LEP data 34 ) and 80 GeV, while letting the gluino mass 
vary between 100 and 200 GeV. In the corresponding range of MSSM parameter 
space, the mass of the LSP (x?) is approximately given by Mg/6. 

For X\ — * W^x®, the 1.6% net branching ratio quoted earlier yields between 
70 and 1 dilepton events (of which half are like-sign) per 25 pb _1 at the Tevatron 
for Mg between 100 and 160 GeV. In contrast, if the X\ decays to iv and/or £u, 
the larger 15% net branching ratio can result in up to 1100 (15) dilepton events 
for Mg = 100 (200) GeV, depending on the decay mode and the various masses. 

To illustrate, we present in Table 2 the dilepton event rates (for an integrated 
luminosity of L = 25 pb _1 ) for various cases which yield ^ 30 events. For compar- 
ison, with our cuts, 12 events are expected from tt production for m t = 140 GeV, 
but all of these have opposite-sign leptons. Thus, for comparable gluino and top 
quark masses, the additional requirement of two isolated like-sign leptons would 



An isolated lepton is defined to be one that is separated by at least 0.3 units in AR = [(A77) 2 + 
(A0) 2 ] 1//2 from any parton (or "merged parton jet" if two or more partons are within 0.7 units in 
AR of each other). 
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reduce the top quark rate to a level far below that from gluinos. However, a top 
quark significantly lighter than the gluino might require that additional cuts be 
made to separate the two signals. 

The event rate depends strongly on the lepton cuts. The lepton spectrum itself 
is quite sensitive to decay modes and decay product masses. For the chargino three- 
body decay to x\^ v Yia - a virtual W , the lepton spectrum depends primarily on 
M~± —M~o. For chargino decays to iv (vi), the spectrum is essentially determined 
by Mj— M~o (M~± — M~). Because of the current experimental limit of Mj> 44 

GeV 34 , most leptons from i decay pass our cuts for the M~o values employed. Thus, 

Al 

even for the relatively large gluino masses of 160 and 180 GeV, the iv decay event 
rates illustrated are large enough to be in possible conflict with observed rates 
at the Tevatron. In contrast, as illustrated in Table 2, the event rate associated 
with X\ ~^ & decays could be very small since small values of M~± — M~ are 
possible, leading to a soft lepton that is unlikely to pass our cuts. However, even 
for Mg = 180 GeV, if the v mass is not close to M~± then the event rate for the 

y Al 

iv mode is large. 

Thus, if very few or no like-sign dilepton events are found after accumulating 
L = 25 pb -1 , then improved limits on the gluino mass (as a function of other 
MSSM parameters) will be attainable. If X\ ~^ W^X? is the dominant decay, a 
modest improvement of Mg > 120 GeV is possible, based on the like-sign dilepton 
search. In contrast, if iv or iv decays of the X\ are dominant and the i spectra are 
not suppressed by a small mass difference, then limits of order Mg ^ 200 GeV will 
be obtained over the large region of parameter space for which BR(g — > qq'xi) is 
substantial. 

Assuming that one has succeeded in isolating gluino candidates, it is important 
to ask if one can estimate the mass of the gluino, and the masses of the decay 
products. If we define a hadronic jet as having Et > 15 GeV and pseudorapidity 
\r)\ < 2.5, the gluino pair events will typically have between 2 and 4 jets. The 
Et spectrum of the jets is completely determined by Mg — M~±. For the extreme 
case of M § = 100 GeV and M~± = 80 GeV, all jets are too soft (E T < 15 GeV) 
for identification. In contrast, the missing Et spectrum is not strongly dependent 
on the mass splittings or chargino decay modes, at least for gluino and chargino 
masses in the ranges considered here, and is centered at about 50 GeV. We have 
performed detailed studies and find that the most useful distributions for mass 
determinations are E T max , the Et of the most energetic lepton and E T ma,x , the 
transverse energy of the jet with largest Et- For any given decay chain, these 
can be used to estimate the mass differences, and the overall event rate can then 
be used to determine the absolute mass scale, provided statistics are adequate. 
However, with L = 25 pb -1 , event rates for W- mediated decays of the X\ at 
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the Tevatron are inadequate. For £u and/or u£ decays, predicted event rates are 
such that ±25 GeV would be achievable for M~ g ^ 160 GeV. But, as already- 
noted, such large event rates are probably inconsistent with current observations 
at the Tevatron. Scenarios consistent with current Tevatron rates would require 
L <; 100 pb _1 to achieve a ±25 GeV estimate of Mg. In this regard, it is important 
to note that knowledge of the relative importance of the iv and vl decays of the 
X\ is, in principle, not required in order to estimate Mg. The Efp max spectrum will 
allow an estimate of how many events per produced gluino have passed the lepton 
cuts, and the overall event rate will then allow an estimate of Mg. The E^ max 
spectrum would then allow us to estimate M~±. Of course, if M~± is known from 

Xi Xi 

other sources (e.g. LEP-II), Mg can be estimated from the E^ m&x spectrum alone, 
without relying on the absolute event rate predictions which will have substantial 
systematic uncertainties. 

The like-sign signature for gg production provides a powerful tool, both for 
discovering evidence for supersymmetry and for estimating the gluino mass. It 
is the Majorana nature of the gluino that yields this striking signature. At the 
Tevatron collider, favorable assumptions concerning gluino cascade decay branch- 
ing fractions, gluino production could yield as many as 1100 dilepton events (after 
significant lepton cuts) in the current 25 pb _1 run, of which half would be like-sign. 
Since large numbers of events are not seen, many new constraints on the masses of 
the g, X\ i X?) £ an d v will be obtained. If the handful of dilepton events currently 
observed at the Tevatron were due to gg production, an integrated luminosity of 
at least L = 100 pb _1 would be required to obtain a ±15% estimate of the gluino 
mass. 

5. Conclusions 

The Tevatron is expected to have accumulated nearly 100 pb _1 of data by 
the end of the 1993-94 run. The data that will be collected by the CDF and DO 
Collaborations will be sensitive to gluino masses up to 200 GeV. The detection of a 
gluino signal may occur through a variety of techniques. In addition to the classic 
missing transverse energy signal, experimenters should search for isolated like- 
sign di-leptons and isolated hard photons in events with jets and missing energy. 
Detection of a positive signal in at least two of these three channels could be critical 
for confirming a gluino interpretation of the detected events. In addition, with a 
signal in more than one channel, one can begin to zero in on specific values of some 
of the MSSM parameters. 

In the absence of any of the signals mentioned above, one would conclude 
that gluinos (and probably squarks as well) are too heavy to be produced at the 
Tevatron in the near future. Although a modest improvement of gluino mass limits 
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could be achieved with the upgraded Fermilab main injector (perhaps, detectable 
gluinos with M~ g < 300 GeV for L = 1000 pb _1 ), one will need the services of a 
supercollider (either LHC or SSC) to definitively establish or rule out the existence 
of the gluino of low-energy supersymmetry. 
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FIGURE CAPTIONS 

1) The branching ratios for g — > qqx® and g — > qq'xf as a function of fi for 
Mg = 120 and 300 GeV and tan/5 = 1.5 and 4. Sections of the curves that are 
not plotted correspond to parameter choices that yield M~± < mz/2. Taken 
from Ref. 15. 

2) BR(x2 ~ * Xi7) as a function of the MSSM parameter //, for tan/3 = 2. The 
three curves correspond to different choices of gluino masses as indicated on 
the figure. 

3) BR(x2 — > Xi7) a s a function of the MSSM parameter for tan/3 = 10. The 
three curves correspond to different choices of gluino masses as indicated on 
the figure. 

4) BR(x2 ~^ Xi7) as a function of tan/3 for two different sets of gluino masses 
and choices of /i. 
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